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Abstract 
O . 
£Sj , Recent measurements in the charm sector at BABAR are reviewed. The scope 

of the presentation includes the observation of two new narrow mesons in the 

-D+7T and Dfir°-f final states as well as the measurement of D°-D° mixing 

parameters by means of two methods: using the doubly-Cabibbo-suppressed 

D° decay to K + ir~ and using the ratios of lifetimes extracted from samples of 

D° mesons decaying to K~ir + , K~K + , and ir~ir + . 
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1 Introduction 

B-factories are a powerful tool for charm physics: the effective e + e~— > cc(bb) 
cross section at the energy of the T(45 l ) is 1.30(1.05)nb. The recent BABAR 
charm results reviewed in this paper concern D°-D° mixing with 2-body hadro- 
nic decays and cs spectroscopy. 

Mixing searches are discussed (Sec. 2) and results from two different 
experimental methods are presented l=» : lifetime difference searches using D° 
Cabibbo-suppressed decays (CSD) to K~K + and it~it + (Sec. 2.1) and wrong- 
sign searches using D° doubly-Cabibbo-suppressed decay (DCSD) to K + tt~ 
(Sec 2.2). The D° Cabibbo-favoured decay (CFD) to K~tt + is needed in both 
cases as reference mode. The former(latter) analysis is based on a 91(57.1) fb~ 
data sample collected on or near the Y(45) resonance with the BABAR detector 
at the PEP-II asymmetric-energy e + e~ storage ring l=». 

The observation of the two narrow states D*j(2317) and D s j(2458) (Sec. 
3) is based on an integrated luminosity of 91fb~ . 

The kinematical selection from continuum is obtained by requiring a mi- 
nimum value of about 2.5 GeV/c for the momentum of a charmed meson in the 
center-of-mass frame. This criterium provides full removal of charmed mesons 
from B decays and a strong reduction of the combinatorial background. En- 
hanced rejection can be obtained by increasing this threshold when needed as 
in the spectroscopic studies of Sec. 3. 

The D°^>h + h~ candidates selected in the mixing analyses presented in 
Sec. 2 are mostly D*-tagged: they are required to be produced by the decay 
D* + ^D°Trf . The slow pion charge provides the flavour tagging; it also allows 
to introduce the quantity 5m = rn(h + h~7r s ) — m(h + h~) that is crucial to 
improve the background rejection. 

Vertexing, tracking and particle identification performances are crucial 
for the analyses presented. The standard kaon identification provided mainly 
by the Cherenkov detector is characterized by an average efficiency of about 
85% and an average pion contamination probability of about 2%. 

2 D°-D° mixing 

The two D° physical states can be represented as |Di,2) = p\D°) ±q\D°) (with 
masses m\^ and widths T\p) where \p\ +|(?| =1 and |D°),|D°) are the flavour 



eigenstates. The mixing parameters x = m\ — rri2/T, y = T± — T2/2T, where T is 
the average width, are predicted by the Standard Model (SM) to be very small 
in the charm sector: \x\, \y\ < 1CP 3 ; however new physics contributions may 
enhance x, whereas final state interactions and SU(3) breaking can enhance 
y "*». Thus charm mixing provides a sensitive probe for physics beyond the 
SM at the level of the current experimental sensitivity (^10 -2 ). 

Since new physics may not conserve CP, it is necessary to consider CP 
violation (CPV) when measuring mixing. CPV effects can be parametrized in 
terms of the following quantities (D°,D° decay into a final state /): 



, <p f = «rg(±^\ ; A f = {f\H D \D°) , A f = (f\H D \D°) (1) 



A value of r m ^ 1 would indicate CPV in mixing; a value of iff ^ would 
indicate CPV in the interference of mixing and decay. 

2.1 Lifetime difference searches 

Mixing would alter the decay time distributions of D° and D mesons that 
decay into CP eigenstates. They can be considered to a good approximation, 
as pure exponential with effective lifetimes "» 

r ± = t° [l + r^ 1 (ycosipfTxsm(p f )] ' (2) 

where r° is the lifetime for the CFD D° -> K-jr+ (and D° -* K+ir~) and r+ 
(r _ ) is the lifetime for the CSD of the D° (D°) into CP-even final states (such 
as K~K + and ir~ir + ). They can be combined to build the parameters 

Y-^-l AY-^A where i (r) = (r + +r~)/2 

Y ~ (r) L ' Ay " <r) At ' Wh6re \A T = (r+ - t-)/(t+ + r") [6) 

Both Y and AY are zero if there is no mixing. Otherwise, in the limit of 
CP conservation in mixing, Y = y cos ipf, AY = xsmcpj; moreover, in the CP 
conserving limit in the interference between mixing and decay, Y = y, AY = 0. 
In these limits Y is not a new physics search parameter but simply a measure 
of the mixing level due to SM. 

The D° mass distribution for the four independent selected samples are 
shown in figd Three of them are D*-tagged: K~ir + measures r°, K~K + 
and 7T _ 7r + both measure (r) and A T . The untagged sample K~ K + , about six 



Table 1: Y and AY results; the first error is statistical, the second systematic. 



Sample 


Y(%) 


AY (%) 


R-R+ 

7T~7T + 

£>*-untagged K~K+ 


1.5 ±0.8 ±0.5 

i.7±i.2±5:§ 

0.2 ± 0.5 1 0.4 


-1.3 ±0.8 ±0.2 
0.3 ±1.1 ±0.2 


Combined 


0.8 ±0.4 + £4 


-0.8 ±0.6 ±0.2 



times wider than the tagged one at the cost of a purity drop from 97% to 68%, 
measures only (r). The mass fits determine the probability of each selected 
D° to belong to a signal event that is used in the proper time fit. Sidebands 
candidates are included in the fit sample to better constrain the background. 
The lifetime for each sample is extracted by an unbinned maximum likelihood 
(UML) fit; the likelihood is a sum of two probability density functions (PDFs), 
one for the signal and one for the background, each differently modelled as 
the convolution of a decay-time distribution and a resolution function. The 
fit results are shown in fig|21 The average (r) and the asymmetry A T are 
calculated from the lifetimes r° and r . Systematics effects on lifetime tend to 
cancel in the lifetime ratio. This analysis is described in detail elsewhere 12*. 

The results, summarized in tab^ are consistent with the absence of mi- 
xing and CPV within the current sensitivity. A comparison with other experi- 
ments is given in tab|2 BABAR provides the most stringent limits on Y and 
the first measurement of a CPV parameter with the lifetime ratio method. 

The Y measurement is statistically dominated (by CSD) and a 1% sensi- 
tivity, for a single tagged sample, is within the reach of the B-factories. 



Table 2: Comparison for Y and CPV results (Belle's Ar is A T of eq.\Q). 



Experiment 


Y{%) 


CPV parameter (%) 


BABAR combined ^ 


0.8±0.4±gj 


Ay = -0.8 ±0.6 ±0.2 


Focus 121 


3.4 ±1.4 ±0.7 


— 


CLEO II 13 


-1.2 ±2.5 ±1.4 


— 


Belle K~K + {prelim.) 19) 


1.2 ±0.7 ±0.4 


A T -- -0.2 ±0.6 ±0.3 


Belle D*-untagged K~K+ ^ 


-0.5±1.0±g;J 


— 



2.2 Wrong sign searches in hadronic decays 

This method implies the study of the time evolution of the DCSD D° — >K + ir~ , 
the so called wrong-sign(WS) decay, looking for a mixing signal followed by 
a CFD: D° -* D° -> if+TT - . Since DCSD and CFD may be characterized 
by different final state interactions, there is an unknown relative strong phase 
between their amplitudes, Sk-tt, not measurable in a mixing analysis. Assuming 
that i,i/Cl and CP is conserved (r m = l, = 0), the time-dependent WS decay 
rate is given by 

„rt 



T ws (t) oc e rt R D + vWR D Tt + -ix' 2 + y' 2 )T 2 t 2 , T = -* (4) 



-{x l2 +y l2 )T 2 t 2 




x 1 = xcosSk-k + ysinSxTi 


y' = —xsiuSk-k 


+ ycos8 K -K 



where R D = gg^g^ , < \ . , ' '''7 " ZZ^'' '■» 



r(D°^ii:+7r-) 

The second and third terms in eq. ^introduce a deviation from a pure exponen- 
tial; the third one is the mixing term whereas the second one is the interference 
term that enhances the mixing signal. The time-integrated WS decay rate is 



Rws = Rd + v'VRd + Rm , Rm = (x u + y' z )/2 (6) 

where Rm is the mixing rate. Additional CPV effects are included by applying 
eqHto D° and D° separately: {R D , x' 2 ,t/'}^{i?^, x' ±2 , y^}. The quantities 
Ad-, Am are related to CPV in the DCSD and mixing respectively, whereas 
CPV in the interference of DCSD and mixing is parametrized by the phase <j>: 



. _ Rp,M Rp,M f x l±: — Cm{x'cos(J> ±y'sin4>) „ 4 / l ± A M 
* ~ R + d,m+ R d,m ; I V* - CM{y'cos4> ± x' Sl n4>) ' ° M V T^Am~ (?) 

Each reconstructed candidate is assigned to one of the four categories 
based on its origin (D° or D°) and its decay (RS or WS). The mixing pa- 
rameters are determined by an extended UML fit to the RS and WS samples 
simultaneously in a four-dimensional space of variables: mass rriKir, mass dif- 
ference 5m, proper time t and its error at- Within a category, the likelihood 
is a sum of PDFs, one for each signal or background component, weighted by 
the number of events for that component; each component's PDF is modelled 
as the convolution of a decay-time distribution and a resolution function. 

The fit is performed in two main steps. Firstly a fit is done in the m^„ — 
5m plane to extract the number of candidates in the signal and background 



Table 3: Wrong-sign results: 95% C.L. intervals including systematic errors. 





95% C.L. range or value (all xlO" a ) 


CPV allowed 


No CPV 


No Mixing 


No CPV nor Mixing 


Rd 


(2.3,5.2) 


(2.4,4.9) 


3.57 ±0.22 ±0.27 


3.59 ±0.20 ±0.27 


A d 


(-2.8,4.9) 




95 ±61 ±83 




x' 2 


< 2.2 


<2.0 






y' 


(-56,39) 


(-27,22) 






Rm 


< 1.6 


<1.3 







components. Its result is shown in fig. [3] for the WS sample; the high level 
background has two peaking components: one for tuk-k, associated to a real D° 
with wrong slow pion, and one for dm, associated to Kir doubly misidentified. 
There are about 440 WS signal decays. Secondly, a simultaneous fit to the 
RS and WS decay time distributions is performed. The larger and clean RS 
sample fixes tq and the resolution model parameters for unmixed decays, while 
the mixing parameters are determined only by the WS sample. The background 
time distributions are determined from ttik-k and 5m sidebands in data. 

A procedure based on toy Monte Carlo experiments is used to set fre- 
quentistic 95% C.L. contours in the fitting parameters (shown in fig. HJ. This 
analysis is described in detail elsewhere LL12J. The results in tab|31are consistent 
with the absence of mixing and CPV. 

A comparison of Rws, x' 2 and y' among experiments is given in tab^J 
Assuming no mixing, the WS time integrated decay rate is just the rate of 
the DCSD. For mixing parameters results, a direct comparison is possible only 
when CP conservation is assumed. BABAR's limits include systematic un- 
certainties and have been obtained allowing x' 2 < 0. The rather different 
techniques used to obtain limits and treatment of the fit parameter x' may 



Table 4: Partial (no CPV) comparison of results of the wrong-sign analyses. 



Experiment 


Rws (96) 


x rz (%) 


y'(%) 


babar 10 > 


0.357 ±0.022 ±0.027 


< 0.200 


(-2.7,2.2) 


Focus (prelim.) ^* 


0.404 ±0.085 ±0.025 


< 0.152 


(-12.4,-0.6) 


CLEO II 121 


0.332 ±g;g|| ±0.040 


< 0.076 


(-5.2,0.2) 


Belle (prelim.) 1^ 


0.372 ± 0.025 ± oloxS 


— 


— 



explain the reason for which, in spite of a data sample six times larger, the 
BABAR upper limit for x' 2 is almost three times larger than CLEO's (two 
times when considering the systematic error ±0.04% quoted by CLEO). More- 
over toy Monte Carlo studies l=^» show that for a given experiment, when 
the preferred sign of y' is positive (negative) as in BABAR(CLF,0) analysis, the 
allowed region becomes larger (smaller). Naively this can be explained by a 
partial cancellation (that takes place with y' < 0) between the interference and 
mixing terms in eq. 0] so that the fit becomes sensitive to small changes of the 
mixing parameters. 

At B-factories the WS analyses can improve x sensitivity towards 1% level. 
Promising techniques have been proposed to measure the relative strong phase. 

3 Two new charmed mesons: L>* 7 (2317) and D s j(2458) 

BABAR has observed two new narrow states in the Dfn 1==* and Dfn°"f l==* 
invariant mass spectra, from inclusive e + e~ interactions, denoted as D*j(2317) + 
and Z? s j (2458) + respectively. Both states have been observed by CLEO l=±» 
and Belle l==» experiments as well. Belle 1=** and recently BABAR 1=1* have 
reported evidence for both states not only in continuum but also in B decays. 

The observation for the D* sJ (231 7) — >D+7r° is given in figs.|3^ and|3jD for 
P*(D+tt°) > 3.5 GeV/c and two different Z?+ decay modes. The observation 
for the D s j(2458) + — ► Dfir°j is provided by the peak in the distribution of 
the quantity Am(D* + Tr°) = m(DfTT°j) — m(Z}+7r°) shown in fig.|SJ:: the peak 
is missing when considering D* sidebands; the residual background would cor- 
respond to Z)*j(2317) mesons associated to random photons selected to be 
consistent with the D* -^D s ^ decay. To disentangle these two possible decays 
and evaluate a possible D*j(2317) feed-up effect, an UML fit using the chan- 
nel likelihood method is applied; this fit is simultaneously performed to the 
-D s 7r°7, D s tt° and D s "f mass spectra of all D s tt°^ combinations. Fit results are 
shown in the three plots at left-side of fig. [fjj No clear indication of decays to 
Z?*j(2317)7 is found (the estimated limit is given in tabJSJ. 

On the other hand a feed-down effect due to D s j(2458) has been investi- 
gated; it would be produced as follows: D sJ (2458) + ^D* + tt° and D*+^D+~/ 
where the D+ir° combination could emulate a D*j(2317) + candidate and the 
residual photon would appear to be unrelated or simply lost. Fig. 01 shows the 
fitted D s ir° mass spectrum in data with the peaking background contribution, 



estimated from Monte Carlo simulation, associated to D s j (2458) decays. 

From previous fits, masses (given in tab|3J, widths and yields are ex- 
tracted; from these yield, once corrected by recostruction efficiency, the rela- 
tive production rate in continuum is calculated (tab|SJ. The observed widths 
( ;$ f MeV) are consistent with the experimental resolution as determined by 
simulation studies; thus the intrinsic width must be below the observed one. 
This analysis is described in detail elsewhere LL» lAl*. 

Further experimental information on the decay modes of D* sJ (23 17) and 
Dsj (2458) (including those produced in B decays) allows one to infer their 
most likely spin-parity assignment lsJJ : + and 1 + respectively. Therefore 
they could be identified with the missing j — 1/2 doublet of P-wave cs states. 
Their masses are significantly lower than predicted by established potential 
models (below DK and D*K thresholds respectively) and enough to close off 
the most natural decay modes. If they are interpreted as ordinary cs states, this 
forces them to decay mainly via isospin- violating transition making their widths 
quite narrow; a possible mechanism is the decay through a virtual r\ followed 
by t] — ir° mixing. The une xpec ted masses and widths have led to a great effort 
on theoretical side as well Issr; either these models need modification or the 
observed state might have a non-qq nature. Among the "exotic" explanations, 
a DK molecule, a four-quark state and a Dn atom have been proposed; among 
"ordinary" explanations there have been attempts to revise or modify potential 
models and to couple chiral perturbation theory with HQET; it has been even 
proposed that these states are a mixture of cs and four-quark states. 

Further experimental information on the radiative decays and di-pion 
emission of the new mesons can definitely shed light on their nature. 



Table 5: Summary of D* sJ (2317) and Z? s j(2458); masses are in MeV/c 2 . 

Ibabar^ m 



Measured quantity 



Belle(cont. 



18) 



CLEO 13 



bir+,K* V K+ 



Df modes 



<j)TT^ 



m(D* sJ {23l7)) 



2317.3 ±0.4 ±0.8 



2317.2 ±0.5 ±0.9 



2318.5 ±1.2 ±1.1 



m(£> aJ (2458)) 



cr-g(A,,j(2458)^D;7r u ) 
B(D s j{2458)^D s j) 



2458.0 ±1.0 ±1.0 



2456.5 ±1.3 ±1.3 



2463.0 ±1.7 ±1.2 



0.25 ±0.04 



0.29 ±0.08 



0.44 ±0.13 



g(D sJ (2458)'->D*7r») 

g(g s j(2458)^D s i J (2317) 7 y 
B(D s j(2458)->P;7r°) 



0.55 ±0.13 ±0.08 



< 0.49 (90% C.L. 



< 0.22(95% C.L.) 



< 0.58(90% C.L.) 
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Figure 1: The reconstructed D° invariant mass distributions (points) and the 
fit result (superimposed curve) for the four samples. The portion of each sample 
assigned by the fit to the background is indicated by the shaded region. 
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Figure 2: The D° proper time distribution within a ±15 MeV/c 2 mass signal 
window (points) and the fit result (superimposed histogram) for the four sam- 
ples. The shaded histogram is the portion assigned by the fit to the background. 
The points beneath each plot show the bin-by-bin difference between data and 
fit divided by the statistical error. 
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Figure 3: Wrong-sign distributions of (a,c) K + ir~ invariant mass and (b,d) 
5m mass difference, for (a,b) signal regions and (c,d) sidebands regions. Data 
are shown as points with the fit contributions superimposed: signal (open his- 
togram), unassociated slow pion background (dark shaded), double misidentifi- 
cation background (black) and combinatorial background (light shaded). 
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Figure 4: Fit results and 95% CL. contours in the (x' 2 ,y) plane. The solid 
point represents the most likely fit result assuming CP conservation and the 
open circle the same but allowing CPV and forcing x' 2 > 0. The dotted (solid) 
line is the statistical contour for CP conservation (CPV allowed); the dashed 
(dash- dotted) line is the corresponding contour when systematics are added. 



450 


^ 






i <o 


"o 400 
>-350 
2 300 


L 


f 


I 


' 


lD 250 
<£ 200 
£ '50 


r 




'''V 




S 100 


r 






^wvviwv, 


50 


= 




I , , , , 


i 



2.1 2.2 



2.3 2.4 2.5 



140 
i 120 
1*100 




2.3 2.4 2.5 



100 



80 



> 

2 60 

>40 



20 





c) 



0.25 0.5 

Am(D* + s 7T°) GeV/c 2 



m(D. 7i°) GeV/c 2 

Figure 5: TTie Dfir mass distribution for the decays (a) Df — » i^ + ii' _ 7r + 
and (b) Df -^K + K~tt + it ; the fit results (curves) for the D*j(2317) + signal 
are superimposed. The very narrow peak near threshold is associated to the 
D* + — > Dfir a decay, (c) The Am^o distribution for D fir j combinations in 
the D^ + signal region (points) and sidebands (shaded histogram) . 
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Figure 6: Fit results from channel likelihood method (curves) overlaid on 
the Dfir°~f mass distribution with (a) no weights and after applying weights 
corresponding to the decays to (b) D* + ir° and (c) D*j(2317) + 7. (d) The 
Dfir mass spectrum (with no additional 7 requirement) . The fit results are 
superimposed; the dashed and lower solid curves show the contributions from 
-D s j(2458) + decays and combinatorial background respectively. 



